Astronomy Reports 46, 366-379 (2002) Russian version in Astronomichekii Zhurnal, 79, 407-421 (2002) Printed 15 May 2009 



Non-Steady State Accretion Disks in X-Ray Novae: 
Outburst Models for Nova Monocerotis 1975 and Nova 
Muscae 1991 



G. V. Lipunova^*, N.I.Shakura^'^ 

^ Sternberg Astronomical Institute, Universitetskiy pr., 13, Moscow 119992, Russia 
"^Max Planck Institut for Astrophysics, Garching, Germany 



Received June 25, 2001; in final form September 13, 2001 



ABSTRACT 

We fit outbursts of two X-ray novae (Nova Monocerotis 1975 = A 0620-00 
and Nova Muscae 1991 — GS 1124-683) using a time-dependent accretion disk 
model. Tlie model is based on a new solution for a diffusion-type equation for 
the non-steady-state accretion and describes the evolution of a viscous a-disk 
in a binary system after the peak of an outburst, when matter in the disk is 
totally ionized. The accretion rate in the disk decreases according to a power 
law. We derive formulas for the accretion rate and effective temperature of the 
disk. The model has three free input parameters: the mass of the central object 
M, the turbulence parameter a, and the normalization parameter 8t. Results 
of the modeling are compared with the observed X-ray and optical B and V 
light curves. The resulting estimates for the turbulence parameter a are similar: 
0.2-0.4 for A 0620-00 and 0.45-0.65 for GS 1124-683, suggesting a similar nature 
for the viscosity in the accretion disks around the compact objects in these 
sources. We also derive the distances to these systems as functions of the masses 
of their compact objects. 
DOI: 10.1134/1.1479424 



1 INTRODUCTION 

Accretion provides an efficient mechanism for energy re- 
lease in stellar systems, making many astrophysical ob- 
jects observable. If the matter captured by the gravitation 
of the central body possesses nonzero angular momentum 
relative to this body, accretion occurs in a disk. This is 
true, for example, in binaries, where the angular momen- 
tum is associated with the orbital rotation of the compo- 
nents. In the course of accretion onto a compact object 
whose radius is comparable to the gravitational radius, 
a substantial fraction of the total energy of the accreted 
matter mc? is released. 

Outbursts reflect one of the most fundamental prop- 
erties of accretion: its non-steady-state character. Cur- 
rently, a number of different models are proposed to ex- 
plain non-steady-state processes in accretion disks and to 
describe the observed source variability. One problem is 
to find an adequate description for the viscosity in the 
accretion disk: viscosity is essential for the accretion, and 
the viscosity characteristics specify the features of time- 
depen dent disk behavior. 

In lLipunova fc Shakural (|2000l ') , a new solution for the 
basic equation of time-dependent disk accretion is found 
and applied to a model of an accretion a-disk around a 



compact star in a close binary. An important property 
of the disk in a binary system is that its outer radius 
is limited. Angular momentum is carried away from the 
outer boundary of the disk due to tidal forces, so that the 
rotation of outer parts of the disk is synchronized with 
the rotation of the secondary. It is assumed that the size 
of the accretion disk specified by the tidal interactions is 
constant over the time interval considered. Another as- 
sumption is that the rate of mass transfer from the sec- 
ondary to the disk is small compared to the accretion rate 
within the disk. 

The last condition is satisfied, for example, in an 
X-ray nova outburst. The accretion rate in the disk 
during the outburst reaches tenths of the Eddington 
rate 10"^ (i^M/Mq) Mo/yr or more, where 9 is the 
accretion efficiency and M is the mass of the cen- 
tral object), while the observed rate of mass trans- 
fer from the companion in quiescent periods is 10~^^- 
10"^ ^ Mg/yr (jTanaka fc Shibazaki|[l99g : ICherepashchukI 
X-ray novae are low-mass binaries contain- 
ing a black hole or a neutron star (see, for exam- 
ple , |ChiiEiih^ llooO) ) . The other component, a low- 
mass dwarf, fills its Roche lobe, so mat ter continuously 
flows into the disk (|Cherepashchukll2000l ). 
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Cur rently, more t han 30 X-ray novae are 
known l|Cherepashchukl I2OO0I) . Most of them have 
hght curves with simila r exponentially decreasing pro- 
files IChen et al.1 Il997l ). During the burst rise, the 
intensity increases by a factor of 10^-10^ over several 
days, whereas the exponential decrease of the light curve 
lasts for several months, with a characteristic time of 
about 30-40 days. 

Two mechanisms for X-ray nova outbursts are de- 
veloped: disk instability and unstable mass transfer from 
the secondary. A final choice between them has not been 
made, and each model faces some problems (see, for ex- 
ample, ICherepashciiu^ (|2000l )). In disk-instability models, 
during the outburst, the central object accretes matter ac- 
cumulated by the disk over decades of the quiescent state. 
This idea is supported by the fact that the mass-transfer 
rate in quiescent periods is comparable to an accretion 
rate corresponding to the outburst energy divided by the 
time between outbursts (T anaka fc S hibazaki 1996). 

In any CEise, an important point here is that we as- 
sume the presence of a standard disk at the time of maxi- 
mum brightn ess of the source, as confir med by spectral 
observations l|Tanaka fc Shibazakil Il996l ) . By " standard 
disk" we mean a multi-color a-disk whose inner radius 
coincides with the last stable orbit around the black hole, 
with the velocity of radial motion of gas being small com- 
pared to othe r characteristic velocit i es in the disk. Then, 
the solution of lLipunova &: Shakural (|2000l ) can be applied 
to the decay of an X-ray nova outburst. The light curve 
calculated using the solution describes observations when 
the contribution of the accretion disk dominates in the 
soft X-ray radiation of the system. This phase is char- 
acterized by a definite spectral state and can be distin- 
guished in the evolution of an X-ray nova. 

Modeling the light curves of an X-ray nova in several 
X-ray and optical bands enables one to derive the basic 
parameter of the disk, the turbulence parameter a (this 
is a new, independent method for determining a in as- 
trophysical disks), as well as the relationship between the 
distance and mass of the compact component. 

Here, we apply our model to Nova Monocerotis 
A 0620-00, which is the brightest nova in X-rays observed 
to the present time, and Nova Muscae GS 1124-683.0 

Currently, the most likely mechanism for turbu- 
lence and angular-momentum transfer in accretion disks 
is thought to be Veli khov-Chandras e khar m agnetic- 
rotat ional instability (|Velikhovl 1 19591 : IChandrasekhan 
Il96ll ). w hich is investigate d in a pplication to accretion 
disks by iBalbus fc HawlevI l|l99lh . Calculations suggest 
that this type of instability corresponds to a ~ 10" 



The parameter a, which was introduced by IShakural 

l|l972h . describes large-scale turbulent motions. The large- 
scale development of M HD turbu le nce h as be en sim- 
ulated , for example, by lArmitag^ (|l998l ) and iHawlevI 
(|2000l ). and it was estimated that a? ~ 10 

On the other hand, the following estimates for a were 
derived by comparing theory and observations: for dwarf 
novae ~ 0.1 during an outburst and ~ 0.02 in a quies- 
cent state in a model with a limit-cycle instability (see, 

^ Observations collected from the literature and 
presented in uniform units can be found at 
[htt p : / /xray. sai . msu. ru/ ~ galj a/xnov/ ( 



for example. ICannizzo et al. (Il988l)l; ~ 10~^ fo r disks in 
galactic nuclei Jsiemiginowska fc Czernvl 1 19891): ^ 1 for 
Sgr A * in an advection-dominated model l|Naravan et al.l 
119951 1: and ~ 0.1 -0.3 in the inner, hot advective 
part of the disk for the X-ray novae GS 1124-683, 
A 062 0-00, and V404 Cyg , basing on spectra in the low 
state iNaravan et al.|[r996l '). 



2 MODEL FOR ACCRETION DISKS IN 
X-RAY NOVAE 

The evolution of a viscous accretion disk is descr ibed by 
the d iffusion-type nonlinear differential equation l|Filipovl 
|1984| 1: 



dF 
'dt 



D 



IF 'dh? 



(1) 



where F 

forces acting between adjacent rings of the disk divided 
by 2-K, Wrtp is the component Wrip of the viscous stress ten- 
sor integrated over the thickness of the disk, h = V GMr 
is the specific angular momentum, and M is the mass of 
the central object. The dimensionless constants m and n 
depend on the type of opacity in the disk. If the opacity is 
determined largely by absorption (free-free and bound- 
free transitions), then m — 3/10 and n = 4/5. The "dif- 
fusion coefficient" D specified by the vertical stru cture of 
the disk relates the surface density E, F , and h (|Filipovl 
ll984l : lLvubarskii fc Shakuralll98'it ): 



{GMf F'^ 



2(1 



(2) 



Relation ([2]) is derived from an analysis of the vertical 
structure of the disk. 

A class of solutions f or Eq. ([l| is derived 
in lLvubarskii fc Shakural |l983) during studies of the evo- 
lution of a torus of matter around the gravitating center 
under the action of viscous forces, which are parameter- 
iz ed by the turbulent viscosity parameter a introduced 
m IShakural(|l972l ). In particular, a solution was obtained 
for the stage when the torus has evolved to an accretion- 
disk configuration, from which matter flows onto the cen- 
tral object. When the accretion rate through the inner 
edge decreases, the outer radius of the disk simultane- 
ously increases — the matter carries angular momentum 
away from the center. In this model, the accretion rate 
decreases with time as a power law. The power-law index 
depends on the type of opacity in the disk. 

An important property of a disk in a binary sys- 
tem is the cutoff of the disk at the outer radius, 
in the region where the angular momentum is car- 
ried a way due to the o rbital motion (see, for ex- 
ample, llchikawa fc Osakil ()l994l )1. Taking into account 
the corresponding boundary conditions, we obtained 
a new solution for Eq. ([T l in a general form for a 
disk with uniform opacity ( Lipunova fc Shakural [2000) . 
We described the vertical structure using calculations 
of lKetsaris fc Shakural (jl998l ) in the framework of the gen - 
erally acce pted a-disk model fehakura fc SunvaevllT973l ) . 
Following iLvubarskii fc Shakura (1987), we considered 
two opacity regimes: with the dominant contribution to 
the opacity made by the Thomson scattering of photons 
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on free electrons, and by the free-free and bound-free 
transitions in the plasma. As a result, we obtained explicit 
expressions describing the time variations of the physical 
parameters of the accretion disk. The solution describes 
the evolution of the accretion disk in a binary during the 
decay of the outburst, while the matter in the disk re- 
mains completely ionized. The accretion rate decreases 
with time according to a power law; however, in this 
case, t he power-law index is la r ger th an that for the solu- 
tion of iLvubarskii fc Shakural (|l987l '): —5/2 compared to 
— 19/16 when Thomson scattering dominates and —10/3 
compared to —5/4 when absorption dominates. The solu- 
tions in the two opacity regimes join smoothly, providing 
a basis for applying a combined solution to describe the 
evolution of a disk with a realistic opacity. Our study of 
disks in stellar binary systems indicates that the second 
opacity regime is realized during times investigatecfl. 

The law for the variation of the accretion rate can 
be easily derived from the condition of mass conserva- 
tion in the disk. Let a Keplerian a-disk in a binary have 
fixed inner and outer radii rin and rout- The mass of the 

"-out 

disk is Md = f 2 7rrEdr. If this mass varies only due 

''in 

to accretion through the inner disk boundary (i.e., the 
accretion rate onto the outer boundary is substantially 
lower), then dMd/dt = Min(t). Let the disk parameters 
F and S be represented as the products F{t) /(^) and 
E(t)(T(5), where /(^) and (t(^) are dimensionless func- 
tions of the radial coordinate ^ = h/ho- The value of 
ho is fixed and equal to the angular momentum at the 
outer disk boundary. Let us express the surface density E 
in the integrand in terms of F using Eq. ([2}. It follows 
from the equation for the angular momentum transfer 
that Min(t) = -2-K Fit) f'{C)/ho (|Lipunova fc Shakural 
I2OOOI ): therefore we obtain 



Table 1. Input model parameters. For parameters denoted 
with "+" we use observed values. 



Md = - ML 



r dr 



= (dMd/dt)^"™ X const, (3) 

where dimensionless y{£^) = f'{^)- We will integrate 
Eq. Q over t. Using the fact that m = 10/3 for the 
Kramers law, we conclude that Ma oc {t + S t)~'^^^ and 
Min (X (t + S f)^^"''"^. Note that, to produce such time de- 
pendence, it is sufficient that the "diffusion coefficient" D 
is constant in time, being a function of the radius. 

The bolometric luminosity L = OMi^c? varies accord- 
ing to the same law as th e accretion rate through the inner 
disk radius. As shown in iLipunova fc Shakural (|2000l ). the 
exponential light curve of an X-ray nova can be then ex- 
plained if the considered spectral band contains the fiux 
integrated over the exponential falloff in the Wien section 
of the disk spectrum rather than a fixed fraction of the 
bolometric luminosity. 

In the model considered, the accretion rate depends 
on the distance from the center. In the central regions, 
the accretion rate is virtually independent of the radius; 
but this is not true for the outer disk. The bolometric fiux 



* Note added to the astro-ph version: we point out that im- 
portant here is the type of opacity in the outer disk, where the 
bulk of the mass is contained. 



(1) 


(2) 


(3) 


M 


Mass of the central object 




Mo 


Mass of the optical component 


+ 


P 


Orbital period of the binary 


+ 


f{Mo) 


Mass function of the optical component 


+ 


a 


Turbulence parameter in the disk 




Nm 


Number of H atoms per cm^ to the source 


+ 




Molecular weight of gas in the disk 0.5 




St 


Inner parameter of the model 





from regions of the disk between rings with radii 0.1 rout 
and Tout is about 6% lower than the values obtained for 
a stationary standard disk model. The optical fiux from 
the time-dependent disk differs from that from a station- 
ary disk by a smaller amount, depending on a size of the 
region producing the optical fiux. This size is specified by 
the mass of the central object and the accretion rate in 
the disk. 

To take into account the effects of general relativity 
in the vicinity of the compact object, we use a modified 
Ne wtonian gravitational pote ntial in the form suggested 
bv lPaczvnskv fc Wiital (Il980l ). For a Schwarzschild black 
hole, this potential is 

^ = , (4) 

r - r-g 

where = 2G M/c? . The accretion efficiency 6 for this 
potential is a factor of « 1.45 smaller than that for a 
Newtonian potential. 



3 MODELING PROCEDURE 

3.1 Derivation of Theoretical Curves 

In t he model for a time-dependent disk in a binary sys- 
tem (jLipunova fc Shakura 2000 '). with the opacity in the 
outer disk defined by the bremsstrahlung absorption, the 
variation of the accretion rate in the disk as a function of 
time is given by the formula 

10/3 



M(h,t) = -27r 



l.22Ay{h/ho) I h, 



ho 



,14/5 



D{t + St) 



(5) 



where t is the time, 5t, a normalizing shift, h, the specific 
angular momentum, ho, the specific angular momentum 
of the matter at the outer radius of the disk, D, the con- 
stant in Eq. ([T]), and 

2/(e) ~ 1.43- 1.61 e^'^ + 0.18^^ (6) 

Table [1] presents the input parameters for the model. 
Three parameters are free: the mass of the compact ob- 
ject M, the turbulence parameter a, and the normalizing 
parameter St. A change of the mass of the optical com- 
ponent Mo affects the disk size, which affects the rate of 
variation of the accretion rate in the disk. The following 
values can be obtained from the input parameters: 

(1) The system semiaxis a, calculated as 



G{M + Mo)P'' 
47r2 



1/3 



(7) 
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assuming that the orbits are circular. 
(2) The system inchnation 



/(Mo)(l+g)2 



1/3 



(8) 



where the mass ratio is g = M /Mo. 

(3) The radius of the last nonintersecting orbit 
around the primary, which depends on the mass ratio 
of the binary components and, in general, does not ex- 
ce ed 0.6 of the Ro che lobe size (values are tabulated 
m iPaczvnskil (Il977l )): this corresponds to the radius of 
the outer boundary of the disk rout (jPaczynskil Il977l : 
llchikawa fc Osakilll994l ). 

(4) The diffusion coefficient D appeared in Eqs. ([1} 
and (O: 



D : 



5.04 X fO'''' a''' (m/O.S)"^/'' {M/Mq) Bf 

[e 



^/i»cm^s-i«/^l 



(9) 



where B( = (H^^ 02 nf n4)~^''"' is a combination of di- 
mensionless parameters specified by the vertical struc- 
ture of the disk, cal culated and tabulated as f unctions 
of the optial depth in iKetsaris fc Shakural 119981 ^ Thus, 
B{ depends on the variable optical depth, which can be 
obtained from the disk characterstics at the radius. We 
found that Bf depends stronger on time than on radius, 
and the time dependence is also not very strong, since 
Bf is the combination of the parameters raised to powers 
much smaller than unity. In the modeling, we adopt Bt 
calculated iteratively at the half-radius of the disk and at 
the middle of the time interval investigated. 

To calculate the effective temperature in the disk, we 
use the formula 



T\r) 



Lo(r)r 

47r(j 



f - 



dr \^ '■^(j) 

r no) 



(10) 



where /(C) = 1.43 ^ - 0.46 C'^''^ + 0.03 C^ uj{r) is the an- 
gular velocity in the disk (which is Keplerian away from 
the compact object), a is the Stefan-Boltzmann constant, 
and Tout is the radius of the outer boundary of the disk. 
In Eq. ([TOj, M is equal to the value M(0,t) defined by 
Eq. ([5]). The central regions of the disk (where r ^ rout 
and the product of the last two factors in Eq. (|10p yields 
approximately unity) produce the largest contribution to 
the X-ray emission; here, the accretion rate is nearly con- 
stant over radius, and the distribution of the effective tem- 
perature essentially coincides with that in a stationary 
disk. We also take into account non-Newtonian nature of 
the metric around the compact object in the central re- 
gions of the disk. For the Paczynski-Wiita potential (|4]), 
one has 



GM 1 



r r — r^ 



(11) 



We assume that the bulk of the optical flux comes from 



^ In iKetsaria fc Shakural lll998h . Table lb should read 
114 instead of Ha ; t h e 5th column should be ignored. 
In iLipunova fc Shakural l l200d) . there is a misprint in D in 
(26) and (31). 



the disk (at distances r ^ rin), while the radiation from 
the "transition layer" at the outer boundary, where the 
momentum is carried away, is significantly less. 

The spectral fiux density detected by an ob- 
server CMi__be__calculatcd using the following for- 
mula (|Bochkarev et al.i.l98& l 



Ly cosi 



exp(— Ty) [erg/cm^ Hz s] , (12) 



where d is the distance to the system, not set a-priori in 
the model, is the optical depth corresponding to the 
absorption toward the source, and Li, is the spectral lu- 
minosity of two sides of the disk according to the formula 



'out 

L^=2 7r J B{iy,T{r)) 2 TTrdr [erg/cUz], 



(13) 



where B{u,T{r)) is the Planck function. 

The model curves Fi (erg/cm^ s) are calculated by 
integrating Fi, over frequency. At X-ray energies, = 
(j{u) Nm , and the cross section for absorption by the hy- 
drogen at oms a(v) is presented in the for m of an analyti- 
cal spline (jMorrison fc McCamm"onlll983l ). The number of 
hydrogen atoms in the line of sight toward the source A^hi 
can be found in the literature, an d can also be ca lculated 
using the approximating formula (|Zombecklll990l ) 

iVni ~ 4.8 X 10^^ S(B - V) atoms/cm^ mag , (14) 

provided the color excess -E(B— V) is known and assuming 
that the main contribution to the absorption is made by 
the hydrogen atoms of the interstellar medium and not 
by those directly related to the source. 

For the optical spectrum, absorption is taken into 
account after integrating F^ over frequency. We use the 
interstellar absorption A\ — 2.5 r/ In 10 instead of the 
optical depth r^; here t is some effective value for the 
optical band denoted by A. 

We calculate the light curves Fi in a chosen time 
interval t G [fi,f2], and we set f = at the peak of the 
outburst. The accretion rate at the inner radius, bolomet- 
ric luminosity of the disk, and fluxes in specifled spectral 
intervals (the X-ray and the optical B and V bands) are 
calculated with a step of At — 1 day. 



3.2 Comparison of the Model and Observed 
Light Curves 

We use the light curves in one X-ray band and in the B 
and V band. 

Spectral observations of X-ray novae can help to de- 
termine in which timporal and spectral intervals the flux 
is dominated by the disk radiation, with much less con- 
tributions from other components. In general, it is ob- 
viously desirable to use data at the softest X-ray ener- 
gies (below ~ 10 keV), since the typical X-ray spectrum 
of an X-ray nova after the outburst peak is a combi- 
nation of the disk spectrum (modeled as emission from 
a multi-colo r black-body disk) and a harder power-law 
component (|Tanaka fc Shibazakilll99^ . For each source, 
we must consider the spectral evolution and exclude time 
intervals when either the contribution from the nondisk 
components (or power-law components) is substantial or 
the disk spectrum differs appreciably from that for a 
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Table 2. Zero-points, central waveleng ths, and effective widths 
for the optical bands llZombeck| |l99(]|) . The corrected value is 
denoted by * ; the value from lZombeda l IlQQOl) is given in paren- 
theses (see explanation in the text after Eq. 1151 1). 



Band 




A, A 


AA„ A 


U 


8.37 


3650 


680 


B 


8.198* (8.18) 


4400 


980 


V 


8.44 


5500 


890 



multi-color disk model. Nondisk components in the emis- 
sion can alter the slope of the decaying light curve, and 
this slope dramatically affects the resulting values of a. 

We reduce observed stellar magnitudes at optical 
wavelengths to fluxes in erg/(cm^s) using the formula 



F, = AA. 10 



-0.4 mi -A^' 



(15) 



where the zero-points and effective bandwidths AAi 
are presented in Ta ble [H We have c o rrecte d A% values 
in accordance with iBochkarev et al.] l|l99ll ). where the 
disk color indices are calculated using the radial temper- 
ature distribution of a standard disk, taking into account 
the actual transmission functions for the optical filters. 
The values calculated with our code, using rectangular 
band-pass for th e optic al filters, coincide with those from 
IBochkarev et al.l (ligQll ) to within 0.01". To achieve this 
agreement, it has been necessary to correct the zero-point 
in the B band. 

Let us consider the disk flux F"^" observed in some 
X-ray spectral band at some time t G [ii,t2]. We can 
derive the distance d to the source for the chosen model 
parameters (see Table[T]) by comparing observed flux with 
theoretical Fi{t), taking into account the absorption to- 
ward the source and the inclination i calculated with ([8|. 
We assume that the disk is coplanar to the binary orbit. 

Further, we calculate model light curves Fi using the 
derived d. We find the color excess E{B - V)'"°'*'=' that 
is required to obtain agreement with the observed optical 
flux in the B band. Each optical fight curve is then cor- 
rected for this _B(B-V)™°'*''' by substructing value Ax/2.5 
from the logarithm of the optical flux. For this, we use the 
values Rb = 4.2, Rv = 3.2, i?u = 5 (|Zombeckl 'l990') in 
the formula relating the color excess and the interstellar 
absorption: 



Ax = Rx E(B - V) . 



(16) 



Observations show that certain portions of X-ray 
light curves of X-ray novae are very close to exponential 
dependences and, thus, can be described by only two pa- 
rameters. We approximate the observed hght curves with 
straight lines log Fi = art + bi, whose slopes are to 
be compared to those of the theoretical curves at some 
point t G [ti, t2] in the middle of the time interval. For 
three spectral bands, we fit fluxes at this t, taking into 
account the observed color excess i?(B — V) and the color 
of the disk B-V. A set of models is selected to provide a 
certain agreement between the X-ray and optical fluxes 
and slopes of the observed and theoretical fight curves. 
The value of is calculated for the selected models. 

A large scatter of the observed disk color B-V, which 
reaches 0.05 magnitudes (see below), can be caused by 



both observational errors and possible short-term fluctu- 
ations of the disk color (due to precession, the superposi- 
tion of the optical flux from the secondary, etc.) that we 
do not examine here. We also do not take into considera- 
tion the contribution of radiation from the central X-ray 
source reradiated from the outer parts of the disk. 

For some of the hght curves, the overall traces are 
rather close to exponential, but the exponential law does 
not agree with the observations by criterion (at the 
0.05 signiflcance level), possibly due to an underestima- 
tion of the observational errors, a superposition of fluctu- 
ations in emission, and model limitations. 



4 MODELING X-RAY NOVA A 0620-00 IN 
THE 3-6 KEV X-RAY AND B AND V 
OPTICAL BANDS 

4.1 Observations 

4.1.1 X-ray Light Curves 

The nova outburst in Monoceros in 1975 (A 0620-00, 
V 616 Mon) was observed in X-rays with the orbit- 
ing observatories Ariel-5, SAS 3, Salut 4, and Vela 
5B jElvis et al.ll 19751: [Poxsev et al.lll97^:lKurt et al.ll 19761 : 
iKaluzienski et al.lll977l : iTsunemi et al.lll989l ). A 0620-00 
was t he first X-ray transient identified with an optical 
burst (iMatilskv et al.lll975l : lBolev et al."l976''). 

W e take the 3-6 keV data of K aluzienski et al] 
l|l977l ), obtained w i th the Ariel-5 All Sky Monitor. Fol- 
lowing IChen et all (|l997D . we assume that the peak of 
the outburst was on Aug 13, 1975, or JD 2442638.3. 
The corresponding data in Crab flux units are taken 
in the electronic form from the HEASARC (High En- 
ergy Ast rophysics Science Archival Research Center) 
database (|Chen et al.lll997t fl. 

We flt the X-ray light curve in units of photons/cm^s. 
Figure [1] illustrates selection of models according to the 
slope of the 3-6 keV X-ray light curve at t = 30 day. The 
regression line constructed using observations in the in- 
terval t G [20, 40] day has a^ = -0.01502 ± 0.0002 and 
br = 1.816 ± 0.007. The data at t ^ W days cannot 
be ascribed with confidence to the exponential section of 
the light curve caused by the disk radiation. The dashed 
curves indicate the boundaries within which we select the 
models. Within these boundaries, the slopes of the lines 
vary in the range (0.9 — 1.1) a^. Value of x^i calculated 
for the observational data and the lines with such slopes, 
divided by the number of degrees of freedom, does not 
exceed 1.3. 



4.1.2 X-ray Spectrum at the Time of the Outburst 

As has often been noted, X-ray novae in ge neral and 
A 062 0-00 in particular (see, for example, iKuulker^ 
(| 19981 )) display softening of the spectrum during the ini- 
tial decay in the light curv e. For A 0620-00 , this was 
pointed out, for example, by ICarpenter et"aL I (119761) for 
3.0-7.6 keV observations and ICitterio et al.l l| 19761 ). for 
3-9 keV. 



ftp: / /legacy. gsfc.nasa.gov/FTP/heasarc/dbase/misc_files/xray_ 



6 G. V. Lipunova and N.I. Shakura 




the range 10^^-10^^ atoms/cm^. For the color excess , we 
adopt value ^(B - V) = 0.35 ± 0.01 (|Wu et al.111983 ). 



t/3 

o 
o 

& 

X 



10 20 30 40 

Time after the peak, days 

Figure 1. Straight-line approximation of the A 0620-00 X-ray 
light curve. The regression line for the data in the interval 
t G [20, 40] days (solid line) and selection boundaries for the 
theoretical curves (t wo dashed lines) are plott ed. The observa- 
tional data are from iKaluzienski et al.l l|l977l V 



Spectral data with high resolution (10 eV at 2 keV 
and 285 eV at 6.7 keV) obtained with the Columbia spec- 
trometer OSO 8 on board the Ariel-5 spacecraft in Oc - 
tober 1975 were analyzed in lLong fc KestenbaumI l| 19781 ). 
The X-ray continuum of A 0620-00 on September 17-18, 
1975 (the 34-35th day after the peak) was best fit by a 
blackbody model with kT ~ 0.5 keV. 

Note that, at the energies of interest to us (> 1 keV), 
the spectrum of a multi-color, blackbody disk can be ap- 
proximated by a Wien spectrum: 

2 



E-" exp{-E/kTm 



(17) 



where Tmax is the maximum temperature in the disk, 
which occurs at the radius ~ 1.58 nn in the case of the 
Paczynski-Wiita gravitational potential Q. In a Newto- 
nian potential approximation, the radius for the maxi- 
mum temperature is approximately 1.36 nn. 

Number of Hydrogen Atoms/cm^ toward 
A 0620-00 and the Color Excess 

The number of hydrogen atoms was estimated 
from the falloff in the soft X-ray spectrum 
due to t he absorption of X- r ays with energies 
< 1 keV llCarpenter et all 1 19761 : iDoxsev et all Il976l : 



iKurt et al.lll976l) : 3 x 1 0^' - 10^^ atoms/cm^ 

In IWu et al.) (|l983l ). the estimate E{B - V) = 0.35 
was derived using UV observations obtained with the 
ANS satellite. Twenty-five stars around A 0620-00 were 
studied, and a specific absorption curve for the region 
was derived. Using (fT4)l . one finds A^hi ~ 1.7 x 10^^ 
atoms/cm'^. The total number of hydrogen atoms in the 
Galaxy in the li ne of s ight toward A 0620-00 is ~ 4 x 
10^^ atoms/cm^ (|Weaver fc Wilhama.1974 ). 

Obviously, we cannot exclude the possibility of sub- 
stantial and variable absorption of the X-ray radiation 
in the source itself. However, the character and origin 
of this variability are unknown. We model the observa- 
tions for Nova Monocerotis for a set of A^hi values in 



4- 1.3 Optical Light Curves 



Optical observa tio ns fro m lLiutvi ll 19761): Shugarovl ( 1976|) ; 



Duerbeck fc Waited 11976) 



van den Berghl ||1976 _ 

Robertson et all l| 19761 ) : iLlovd et all (|l977l ) are used. We 
construct linear regression fits log F — a^t + to the 
B and V optical data in logarithmic flux units for times 
t £ [0, 47] days using the weighted least-squares method. 
This yields = -0.0079 ± 0.0002, br = -9.675 ± 0.005 
for the B band and = -0.0071 ± 0.0002, 

br = -9.885±0.006 for the V band (see Section[321). The 
reduced values for these fits are roughly 12 and 43 for 
the B and V bands, with 102 and 89 degrees of freedom, 
respectively. This suggests that the adopted errors for 
the optical observations might be underestimated (for 
example, systematic deviations have not been taken 
into account) and/or that our assumed exponential 
(quasiexponential) decrease of the optical fiux does not 
yield a complete description of the observed curves, due 
to various fluctuations and variations of the optical flux 
superimposed on the overall trace of the light curve. 
Nevertheless, for modeling we assume that the basic 
trend of the light curves can be fit by a quasi-exponential 
decay. 

In model fitting, we adopt the color of the disk at 
t = 30 days B-V= 0.24 ± 0.03™, derived from the above 
observational data. 



4.2 Results for A 0620-00 

We compare the theoretical and observed curves for t G 
[20, 40] days after the outburst peak. Table |3] summa- 
rizes the parameters attempted. The number of hydro- 
gen atoms along the line of sight toward A 0620-00 does 
not appreciably affect the results, since the absorption at 
3-6 keV is small. 

Figure [2] presents the results of our modeling for 
the parameters from Table [3] We can see that a lies in 
th e range 0.225-0.375 ( f or s lightly different masses Mo). 
In lLipunova fc Shakur^ lioOlif ). we adopted Mo = 0.5 Mq 
and _B(B — V) = 0.39 and used a broader range of slopes 
than in Fig. [T] so a broader interval for a was obtained 
for A 0620-00 in that study, from 0.3 to 0.5. 

Figure |3] shows the resulting relationship between 
the distance to A 0620-00 and the mass of the black 
hole. The distance to A 0620-00 has been estimated t o 
be from 0.5 t o 1.2 kpc (for example, in lOkd (Il977l): 
Ivan den Berghl (119761): see al so reviews by IChen et al.l 
(|l997l ): ICherepash chuk' (2000')). Figure [3] indicates that 
for d ~ 0.9 kpc {Okc 1972,) the mass of the black hole 
in A 0620-00 is ~ 9 Mm. This mass is consistent with 



previous estimates (see ICherepashchukI l|2000l ) and ref- 
erences therein). Figure |3] presents an example of the 
modeled light curves, with i = 47°, d = 0.66 kpc, the 
bolometric luminosity Lhoi{t = 0) = 0.25 1/Edd , and 
rmax(t = 35) = 0.45 keV (cf. Section gXlJ . The reduced 
for the X-ray light curve for t G [20,40] days is 1.17. 
We are not able to obtain a satisfactory fit to the 
slope of the optical light curves. In principle, steeper op- 
tical light curves can be obtained by taking into account 
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Figure 2. Modeled a and color of the A 0620-00 disk. Results 
are given for the values from Table 3. The mass of the optical 
component varies from top to bottom: 0.3, 0.5, and 0.7 Mq. 
The filled circles satisfy the disk color condition B — V = 0.24± 
0.03 mag for t = 30 days. 
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Table 3. Input parameters for A 0620-00 models. The mass 
of the optical component Mo = 0.5 M(7) , its mass f u nction , 
and the binary perio d are taken from IChen et al ] l|l993); 
ICherepashchukI 1 I2OOOI ). 



Parameter 


Tested values 


M 


5-25 A/0 


Mo 


0.3, 0.5, 0.7 Mq 


P 


0.322 day 


fiMo) 


2.7 Mq 


a 


0.1 - 1 


Nm 


3 X 10^1 -10^2 atoms/cm^ 




0.5 


St 


50 - 250 days 



Ef 10 




20 40 

Time after peak, days 

Figure 4. Example of modeling light curves of A 0620-00 in 
the 3-6 keV, B, and V band. The model parameters are M = 
7Mq, Mo = 0.5Mq, a = 0.3, St = 168 days, and Nui = 
3 X 10^-*^ atoms/cm^. 



standard a-disk model, irradiation is insignificant due 
to the smal l half-thickness of the disk, which is typically 
about 0.03 l|Lipunova fc Shakuralfioool ). Probably, a con- 
tribution to the optical flux from a disk, which is warped, 
should be taken into account; a further study of the gen- 
eration of optical radiation in a time-dependent disk is 
necessary. 



Figure 3. The distance-black hole mass dependence for 
A 0620-00 for the model parameters from Table 3. 



irradiation of a thick or twisted disk. The outer parts of 
the disk intercept some of the X-ray flux from the cen- 
tral regions, causing the effective temperature of the outer 
disk, and accordingly its flux, to increase. The intrinsic 
and reprocessed flux depend on the accretion rate in dif- 
ferent ways: the reradiated optical flux decr eases more 
rapidl y, steepening the optical light curves. lEsin et al] 
(|2000l ') have modeled the outburst of A 0620-00 taking 
into account irradiation of the disk and assuming a rel- 
ative half-thickness for the disk of 0.12 (which is appre- 
ciably higher compared to the standard model). In the 



5 MODELING X-RAY NOVA GS 1124-683 IN 
THE 1.2-37.2 KEV X-RAY AND B AND V 
OPTICAL BANDS 

5.1 Observations 

5.1.1 X-ray Light Curves 

The outburst of Nova Muscae (GS 1124-683, GU Mus) 
was discovered independently by WATCH/GRAN AT and 
ASM/ GINGA (All-Sky X-Ray Monitor) on January 9, 
1991 (|Lund et al.lll99ll : [Sunvaev, 19911). The assoc iated 
optical outburst was also detected I Lund et al.lll99ll '). For 
the modeling, we use 1.2-37.2 ke V data obtained wit h 
the GINGA Large Area Counters (lEbisawa et al.lll994 ). 
The data i n erg/(cm^s) are provided by the HEASARC 
database (|Chen et al.1 Il997l . see Footnote |3J. Follow- 



8 G. V. Lipunova and N.I. Shakura 



ing lChen et alj |l993), we take the peak of the outburst 
to be on January 15, 1991, or JD 2448272.7862. 

The weighted least-squares regression line for data 
in the interval t G [35, 61] days yields ar = —0.0134 ± 
0.0001 and = -6.511 ± 0.005. However, the calculated 
reduced is very high, since the observational errors are 
small and the data show an appreciable scatter around 
the general trend of the model light curve. To select the 
models, we use interval (0.98—1.02) for the the slopes of 
lines tangent to the theoretical light curves at t = 48 day. 

5.1.2 X-ray Spectrum at the Time of the Outburst 

As shown by lEbisawa et al.l (|l994l ). after the peak the 
spectrum of GS 1124-683 sof t ened as the luminosity 
decreased. In iKit amoto et al.' (^199 ^: |Mivamoto et al.l 



Table 4. Input parameters for GS 1124-683 models. The 
mass of the opti c al co mponent Mo = 0.8 Mq is adopted 
from IChen et al.l lll997l ). and the binary period and mass 
funct ion of the optic al component 3.01 ± O.15M0 are 
from lOrosz et al.l l ll996l) . 



l|l993h : lEbisawa et al. (1994); Grcin er et al.1 (Il994l '). the 



observed X-ray spectrum was approximated by a model 
with two components: a blackbody, multi-c olor disk and a 
harde r power-law component. Figure 15 in lEbisawa et al.) 
l|l994ll indicates that, in the time interval of interest. Tin 
for the spectral approximation (E bisawa c t al. 1994) is 
roughly 0.7 keV. 

It was suggested bv IKitamoto et aP l|l992l ) that 59% 
of the flux on January 15 (near the outburst peak) was 
blackbody disk radiation, while the remainder was con- 
tributed by a power-law component. During the follow- 
ing 25-30 days, the power-law component decayed more 
rapidly than the disk component. ROSAT o bservations on 
Janu ary 25 (the 10th day after the peak) (|Greiner et al.l 
Il994h suggested that the 0.3-20 keV flux at that epoch 
was completely produced by the disk component. Using 
the approximat ion for the observed 1. 2-37.2 keV X-ray 
spectrum from iMivamoto et al.| (Il993l ) and the derived 
fluxes of the spectral components, we conclude that the 
1.2-37.2 keV flux during t G [35,61] days after the out- 
burst was determined by the disk radiation and can there- 
fore be used in our modeling, since the contribution from 
nondisk components is apparently negligible. 

Number of Hydrogen Atoms/cm^ toward 
G S 1124-683 a nd th e Color Excess 
In ICheng et~ai1 (|l992l ). _B(B - V) ~ 0.29 was estimated 
from HST observa tions of the interstellar abs o rption 
profile at 2200 A. IShrader fc Gonzalez-Riestral (|l993l ) 
found _E(B — V) = 0.3 ± 0. 05 u sing a similar tech- 
nique. In Idella Valle (Il99ll ). the same value 
- V) = 0.30 ± 0.10 was derived from intersteUar 
Na D lines. Using (fT4)) . we arrive at Ai^Hi ~ 1.4 x 10^^ 

ato ms /cm^. 

iGreiner et ai] l|l994h obtained iVni « 2.2 x 
10^^ atoms/cm^ by modeling the combined ROSAT 0.3- 
4.2 keV and GINGA 1.2-37.2 keV data for January 24- 
25 using a composite spectrum with a blackbody, multi- 
color disk and a powerlaw component. For various multi- 
color disk models, they obtained values from 1.7 x 10^^ to 
2.5 X 10^^ atoms/cm^ 



5.1.3 Optical Light Curves 

We use the observational data from iKing et al.1 l|l996l ): 
IdeUa VaUe et al.1 (|l998D and derive weighted least-squares 
regression lines for the optical B and V observations in the 



Parameter 


Tested values 


M 


5 - 25 Mq 


Mo 


0.8, 0.9 


P 


0.433 days 


}{Mo) 


3Mq 


a 


0.1 - 1 


Nm 


(1.4 - 2.5) X 10^1 atoms/cm^ 




0.5 


St 


50 - 250 days 



logarithmic flux units in the time interval t G [12, 61] days. 
This yields a, = -0.0057± 0.0006, b, = -10.79±0.03 for 
the B band, and = -0.0052±0.0006, br = -10.98±0.03 
for the V band. The corresponding reduced values are 
2.8 and 7.3 for B and V, respectively, with 17 degrees 
of freedom. This leads us to a conclusion similar to that 
for A 0620-00 optical light curves (see Section gXS} . We 
again assume that the overall trend of the light curves can 
be described as quasi-exponential. 

In model fitting, we use the observed color of the disk 
for t = 48 days: B-V= 0.27 ± 0.07™. 



5.2 Results for GS 1124-683 

We compare theoretical and observed curves in the in- 
terval t G [35,61] days after the peak of the outburst. 
Figure [5] presents the results of our modeling for the pa- 
rameters from Table |4] Comparing the lower and upper 
left graphs, we can see a slight dependence of the results 
on the number of hydrogen atoms toward GS 1124-683. 
The resulting a values for the GS 1124-683 disk lie in the 
range 0.475-0.625 (for small variations in Mo and A^hi). 

Figure |6] displays the dependence of the distance to 
GS 1124-683 on the black hole mass. Estimates of the dis- 
tance to GS 112 4-683 in the literatur e range from 1 to 8 
kpc. For 3 kpc (|Cherepashchukl |2000| ) . the implied mass 
of the black hole in GS 1124-683 is ~ 8 (Fig. |6l). The 
values for the black-hole mass obtained by us are c onsis- 
tent with the observations (see lCherepashchukI l|200g) and 
references therein). 

Figure [7] presents model light curves for a particular 
choice of parameters. In the corresponding model i = 68°, 
d = 2.6 kpc, the bolometric luminosity I/boi(i = 0) = 
0.47 LEdd, and rmax(t = 48) = 0.44 keV. The model 
satisfactorily reproduces both the X-ray and optical-light 
curves of GS 1124-683. 



6 CONCLUSION 

We have modeled outbursts of two X-ray novae, 
A 0620-00 and GS 1124-683, using the model of a 
ti me-dependent g-disk in a binary system developed 
in lLipunova fc Shakural (|2000l '). 

The turbulence parameter a is estimated as 0.2-0.4 
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Figure 5. Modeled a and color of the disk of GS 1124-683 
for the parameters in Table 4. The mass of the optical com- 
ponent and the number of hydrogen atoms along the line of 

sight toward GS 1124-683 vary from top to bottom: Mo = 
O.SMq, A^hi = 1-4 X lO'^"'" atoms/cm'^ (two upper graphs), 
Mo = 0.9Mq , A^Hi = 1-4 X lO^^atoms/cm^ (two middle 
graphs), and Mo = O.SM©, Nyii = 2.5 X lO^iatoms/cm^ (two 
lower graphs). The filled circles satisfy the disk color condition 
B -V = 0.27 ± 0.07 mag for t = 48 days. 
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Figure 6. The distance-black hole mass relation for 
GS 1124-683, modeled for the parameters from Table 4. 
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for A 0620-00 and 0.45-0.65 for GS 1124-683. The esti- 
mates of a are close to each other, indicating the common 
nature for viscosity in the accretion disks around com- 
pact objects, at least in these two sources. The value of a 
(< 1) points out that the gas in the disks is appreciably 
turbulent. We have also obtained relations between the 
distances to the systems and the masses of the compact 
objects. 

Thus, it is possible for the first time to model ex- 
ponentially decaying light curves of X-ray novae using a 
model with a thin accretion disk with constant a and es- 
timate value of Q. 

Note that the optical B and V fluxes from 
GS 1124-683 can be explained as the disk radiation gen- 
erated locally due to viscous heating. It is likely that in 
the case of A 0620-00, we must also take into account 
the contribution of reradiation by a warped disk. We will 
consider this problem in a subsequent paper. 
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